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Abstract 

This paper examines some probabilistic properties of the class of periodic GARCH processes (PGARCH) 
which feature periodicity in conditional heteroskedasticity. In these models, the parameters are al- 
lowed to switch between different regimes, so that their structure shares many properties with periodic 
ARM A process {PARMA). We examine the strict and second order periodic stationarities, the ex- 
istence of higher-order moments, the covariance structure, the geometric ergodicity and /3— mixing of 
^ ' the PGARCH{p,q) process under general and tractable assumptions. Some examples are proposed to 

■ illustrate the various concepts. 
OS ■ Keywords. Periodic GARGH Processes; Periodic Stationarity; Geometric Ergodicity; Higher- 
. Order Moments. 
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^1 Introduction 

'j_jConsider a periodic GARCH{pi, ...,Ps, qi, Qs) process (xj)^^^ with period s > and ordersp = (pi, 
■ —and q = {qi,...,qs), defined on some probability space {^l,A,P) with the non-linear periodic difference 
equation: 
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where (et)tez ^ sequence of independent and identically distributed (i.i.d.) random variables defined on 
the same probability space (O, A, P) such that E {st} = E {e^} = 0, and E {e^} = 1 (these conditions are 
obviously satisfied if (et)fg2 is Gaussian). 

In the difference equation Xst+v refers to xt during the v — th 'season', 1 < t; < s of period 

t, ao(f ), ai(f ), oip^{v) and /3;^(f ),..., /J^^ (w) are the model coefficients at season v such that for all 
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V = ao{v) > 0, ai{v) > 0, i = l,...,pv, and (3j{v) > 0, j = l,...,qv Moreover, we assume that Sk 

is independent of xt for k > t. We use the periodic notations (xst+v) , i^st+v) , (hst+v) , ioii{v),0 < i < Pv), 
and (/3j(f ), 1 < i < qv) to emphasize the periodicity in the model. There is no foss of generahty in taking 
Py and Qv to be constant in v. If p^ or change with v, one can set p = max Pv, q = max q^ and take 

1<V<S 1<V<S 

Oik{v) = for Pv < k < p and (3j.{v) = for g^, < A; < g, so in the sequel, we shall consider the periodic 
GARCH with constant orders p and q. Since Bollerslev and Ghysels (1996), this type of non-linear models 
has become an appealing tool for investigating both volatility and distinct seasonal patterns, and has been 
applied in various disciplines such as finance and monetary economics (see e.g. Bollerslev and Ghysels, 1996 
and Pranses and Paap, 2000). 

When we consider a periodic model as a data generating process, it is important to find conditions 
ensuring the (periodic) stationarity, ergodicity and the existence of higher moments for further statistical 
analysis. Various probabilistic properties of standard GARCH models have been studied extensively by 
many authors (see e.g., Chen and An, 1998, Bougerol and Picard, 1992a, 19926 and Carrasco and Chen, 
2002 and the references therein). In the present paper, we focus on studying the fundamental probabilistic 
properties of the PGARGH process (xf)^^^ generated by (II. ip so, in Section 2, we present a vectorial 
representation from which we derive some sufficient conditions for the strict stationarity. In Sections 3 
and 4, necessary and sufficient conditions for the second order stationarity and the existence of higher 
order moments are given. Section 5 is devoted to covariance structure. In Section 6 we provide conditions 
under which strictly stationary solutions are exponential /3-mixing with finite higher order moments. We 
conclude in Section 7. 

Some notations are used throughout the paper: /(fc) denotes the identity matrix of order k and 0{kxi) 
denotes the matrix of order k x / whose elements are zeroes, for simplicity we set 0(^j^^ := 0(;jxfc)i Pi^) 
refers to the spectral radius of a square matrix A, i.e., the maximum eigenvalue of a matrix A in absolute 
value, Vec{A) is the usual column-stacking vector of the matrix A, \\.\\ refers to the standard (Euclidean) 
norm in or the uniform induced norm in the space Ai{n) of n x n matrices, denotes the Kronecker 
product of matrices, and A^"^ = A ® A ... ^ A (m— times), for any integer m > 1. For any p > 1, 
JJ' = (f2^ A., P) denotes the Hilbert space of random variables X defined on the probability space 
{0,,A,P) such that \\X\\p = {E \X\^}^^'^ < +oo. We also use the following property of matrix operation, 
Vec{ABC') = {C ® A)Vec{B), where ' is the matrix transpose. 
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2 The Markovian representation and strict stationarity 



Let (xt)^g2 be a process conforming to the model (jl.ip . Setting t/st+v = x'^t+v Vst+v = ^'st+v^ obtain 
from (II. ip the following representation 



(2.1) 



i=l 



Equation (j2.ip is intractable when we want to examine the probabilistic structure of this representa- 
tion. Instead, we will work with the corresponding state-space representation. Let d = p + q and define 
= {Vst+i,--,Vst+sy, y^t+y = {yst+v,--,yst+v-p+i,hst+v,--,hst+v-q+iy and (^7?^^ = {ao{v)r]si-^^,0, ... 
,0, ao(^^)!0) •••jO)' as vectors in M*, and respectively, and set 
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are pxp and pxq matrix valued polynomial functions of r] and where 
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qxp \ / qxq 

this means that Ay{x), By(x) and B^y (x) have entries and coordinates, respectively, which are polynomial 
functions of the coordinates of x. Using the notations above. Equation ()2.ip can now be written as 



yst+v = ^v(.!lt)yst+v^l+^-^!lt^ 



(2.2) 



with Ust+v = H'y^^^^ where H' = (1, 0, 0)^^^^. Equation (j2.2p . is the same as the defining equation for 
multivariate generalized periodic autoregressive process introduced recently by Franses and Paap (2000). 
However, since Gladychev (1961), with periodic time-varying coefficients, it is possible to embed seasons 

is a generalized 

(2.3) 



into a multivariate stationary process. More precisely {Y_t)i^i where y_t= [V 
autoregressive process, i.e.. 
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where A{i] ) and B_{r] ) are defined by blocks as 



/ O 



(d) 



0(d) Mvt) \ 
0{d) M!it)Mvt) 



s-l 



o. 



id) 



v=0 



J 



dsxds 



k=i I i.=o -* ' 



dsxl 



where, as usual, empty products are set equal to /(^^^ . 

In this section, we are interested in strictly stationary and causal solutions for The results of 

this section are based on theorems proved by Bougerol and Picard (1992o) for generalized autoregressive 



representation. Since (r]^) is an i.i.d. process, {A{riJ,B_(r] 



'--it' 



is a strictly stationary and ergodic 



process and since jlog"*" A^rj^) | < | A{ri^) | and -Bllog"*" B_{r]^) | < | B_{rj^) |, then both 
E |log^ ^(^g) I ^iid E |log^ B_{r]^) | are finite where for any 2; > 0, log^ x = max(logx,0). 

Theorem 1 [Strict stationary solution] Equation l\2.3^ has a unique strictly stationary and ergodic so- 
lution if and only if the top Lyapunov exponent {A) associated with the sequence matrices (A{ri^ 



-ij (A) := inf < 
^ ^ t>o ' 



^-log 

t ^ 



t-i 



j=0 



lim { \ log 



t-1 



i=o 



(2.4) 



is strictly negative. The unique stationary solution is causal, ergodic and given by 

00 I A;— 1 



k=l \j=0 



(2.5) 



where the above series converges almost surely (a.s). 

Proof. The theorem is a multidimensional extension of Theorem 1.3 by Bougerol and Picard (1992a). 

Remark 1 Since (A) is independent of the norm, thus we can work with some norms that make it 
rather straightforward to show [A) < 0. However, sufficient conditions which ensure 7^ [A) < are 



1. i?{log||A(7?Q)||} <0. 



2. E 



t-i 



i=o -* ^ 



< 1 for some r > and t'>l, in which case E{\\Y_^\\'^} < +00. 
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Now, a simple computation shows that 



( O 

o 



j=0 



(d) 
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id) 
id) 
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id) 
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Oid) nrnV.-.(!z,_.)| 



Therefore, because the top-Lyapunov exponent is independent of the norm, by choosing a multiphcative 
norm it is straightforward to show that 



7l (^) <7lW :=m.{lE^log 



4 = 1 



Lt-i 



where $ ^ry^^ •= 11 We have thus shown our first result which gives a sufficient condition for 

strict stationarity. 



Theorem 2 Suppose that 7i(<l>) < 0. Then for allt^X, the series 

oo ( k-1 ^ 

k=0 {j=0 ) 

converges a.s.and the process {Y_t)i^j^ defined by l\2.5\i is the unique strict stationary, causal and ergodic 
solution of (fOj) . 



Example 1 In PG ARCH {1,1) , writing (j)^ yrj^j = (^r]g^_^_y,ly {ai{v), Pi{v)) , we obtain 
s-i rs-i ^ 



v=0 



Hence 



log 



11'^ fe-. 



1=1 



s-l t 



i=l 



= YjYj^^^ [Vsit~i)+vai{v + 1) + Pi{v + 1) j + log 

11=1 i=l 

By the law of large numbers, a sufficient condition which ensure 7i(A) < is 

s 



nU(*-)+-i) («i(i),/3i(i)) 



v=l 



which reduces to the classical condition when s = 1. It is worth noting that the existence of explosive 
regimes (i.e., ii^jlog {Vst+v'^ii'^) +/^i(^))} > ^) ^'^^s not preclude (periodic) strict stationarity. 
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The top-Lyapunov exponent criterion seems difficult to obtain explicitly, however a potential 

method to verify whether or not < is via a Monte-Carlo simulation using Equation (|2.3p . This 

fact heavily limits the interests of the criterion in statistical applications. Indeed, the solution need to 
have some moments to make an estimation theory possible and Condition (|2.4p does not guarantee the 
existence of such moments. Therefore, we have to search for conditions ensuring the existence of moments 
for the stationary solution, for which, the top-Lyapunov exponent will be automatically negative (see 
Remark 



3 Necessary and sufficient second-order stationarity conditions 

In this section we examine the necessary and sufficient conditions ensuring the existence of a unique causal 
periodically correlated (PC) solution to (II. ip . This is equivalent to examining conditions of the existence 
of solutions in of the process (Y_^)^^^ defined in equation (|2.3p . Let A = E |^(t/^)| and B_ = E |p(t/^)| 
be the expectations of the random matrix A{r]^ and the random vector B_{rj^), then A and S are defined 
element-wise as 



o. 
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id] 
id] 
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id] 
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Oid] n (t^s-v 

v=0 



.B 
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2B.1 + B_2 



s ( s—k—1 

E n ' 

k=l I v=0 



where 



A,{1) B,{1) 
A^ B^ 



■>B^ 



dxd 



( 



\ 



O, 



(p-l)xl 

qo(w) 

V 0(Q-1)X1 / 



dxl 



with 1 = (1,1,...,!)' E W. To verify that the process {Y_f)^f-^ defined by (j2.5p is well-defined in L-*^, it 

Y[ ^{v^ •) f — fc) converge to zero in L"*^ (endowed 



is sufficient to show that the coefficients 



7=0 



with any matrix norm) at an exponential rate, as k ^ 00. 



Lemma 1 Assume that 



^s-l 



< 1 



(3.1) 



\v=0 



00 I fc-1 I 

then the series ^ \ Yl ^i^t 1 —^^t ^^''^'^^''"9^^ ^"S- Furthermore, the process (Hj^g^ defined by (j^.jp 
fc=o [j=o ) 

is stationary in and satisfying <\2.3^ . 
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Proof. First, we notice that A{ri^) and B_{r]^) are sequences of independent, non-negative random ma- 
trices and vectors respectively, and ^(^/^ •) is independent of B_(rj^ .) for all i ^ j. Therefore, we have 



k-l 



E ^ Yl ^(^t j) \ —(''It fc-' ~ A^B_. It can be shown that the characteristic polynomial of A is det(/(rf5) 



s-l 



XA) = det(/(,) -An <Ps-v), hence p (A) = p MQ </'s-. • If P Tl ^s~v < 1> then ^ A'' < +00, and 

v=0 \v=0 / \v=0 J k=\ 

00 f ffc-l 11 00 ffc-l 1 

X] S S n ^(^(_ ) r ^(^(_^) \ < +00 and thus s Jl ) r ^(^^.j^.) converges a.s.. This further 



fc=o I li=o 



fe=0 I j=o 



k-i 



implies that < Yl ^iVt •) ( ^iVt i.) converges a.s. to the zero matrix as A; ^ cxd. It is obvious that the 
[j=o ^ J 

process (Y_-i.)^^^ defined in p. 51) is stationary in L^. Furthermore, we obtain 



Yt = B{Vt)+MvM 



00 [ k—l 

mt-^)+E{ll^(!Lt^,)\mt-k^ 

k=2 [j=0 

{00 ( k-l 
mt-i)+E{ll^(!Lt-,J\mt-k-i^ 
k=l [j=0 



Lemma 2 // {\2.3^ admits a stationary solution in L^, then p { Yi 4's-v ) ^ I- Moreover, the stationary 

\v=Q 

solution of (l^.gjl is unique, causal and ergodic. 



Proof. From (|2.3|) . we obtain, by recursion, for any n > 1 

n ( k-l 1 

= mt) + E n ^(!Z*-P m-k) + { n ^(i?..,) ^ ^^—i- 

k=l \j=0 I I i=0 



Since, all A{r]^), B_{i]^) and are non-negative, then 



E{Y,}>^E< 



k=0 



j=0 



■t-k' 



^A'^B. 



k=0 



(3.2) 
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This implies that ^ A^B^ < +00. Hence, hm A'^B^ = 0. Since 



k=0 



k—*oo 



( 



A' 



's-l 



k-l 



0(d) ... n 



o 



(d) 



Ki=0 J 

0(d) 4>2(t^i < n ^^-i 

U=0 



k-l 



0(d) ... 



's-l 

n 

.1=0 



's-l 



then to show that hm A^ = 0, it is sufficient to show that hm < TT (t>s-i \ =0 which implies that 



fc— >oo 



.i=0 



-1 



n 



P I 11 VS-* 

\i=0 



p (A) < 1. To do this we shall prove that 



Cj = 0, i = 1, d 



(3.3) 



where (ei,)]^<^<^ is the canonical basis of M'^, i.e., e„ = ((5i,t,, 52,1,, (^d,?;)' where iJjj = 1 if i = j and other- 

( \ ^ ( 

wise. Because the dxdj—th block of vector A''B_ is (A'^B) ■ = S JJ 



fc-i 



's-i-l 



En 



j = 1, s we deduce that for j = s, 

rs 1 
hm <^ TT 



D=0 J Z=l U=0 J v=0 



n '^-w^f 



's-i-l 

n 



■^^ = 0, l = l,...,s. 



(3.4) 



.1=0 ^ I v=0 J 

Since = 00(^^)61 + ao(f )ep+i for v = 1, s and ao(f ) > 0, we obtain from (j3.4p for / = 1, s 



's-l 



's-l 



's-l-l 



U=0 J I v=0 ) U=0 J I i)=0 J 

and we use the relationships (p^ei = ai{v) (ei + ep+i) + ej+i for i = 1, and (/'t,ep+j = /3j(f ) (ei + Cp+i) + 
ep+j+i for i = 1, g — 1 we obtain 



's-l 



's-«-l 







lim <^ TT 

u=o 

rs-l 

lim <^ TT 



n 

. ^=0 

n 



ei = lim <^ TT (p^^i \ Cs+i^i, I = 1, s 



.i=0 

6_|ep+i = £rn 
v=0 J [. 1=0 



Cp+s+l-Z, I — 1, S 



and hence p.3p holds for i = 1, -.-s and i = p + 1, ...p + s. On the other hand, from ()3.4p we have 

{s-i ^ p-i ^ ^=-1 
i=0 ) [i=0 ) 

rs-i -\ ^ (s-i \ '^-1 

fi — ^'Oo I I K — *'Oo I I 

This concludes the proof of the lemma. To prove the uniqueness, let {Z_-f)^^-^ be another stationary process 
conforming to (|2.3p . Then {Zj)^^^ satisfies an equation similar to (j3.2p . By setting = ~ Y-ti 

m 

we obtain, 'im > I : = H ^{v, )W.t^m-i- Defining W_ = E{W_t], we have W_ = A'^+'^W, since 
A := p{A) < 1 we conclude that W_ = (with probability 1). Hence the uniqueness follows. 

Theorem 3 A necessary and sufficient condition for the existence of unique stationary solution in of 
equation l\2.3t} is that (3.1) holds. Moreover, this stationary solution is causal and ergodic. 

Proof. The proof follows from Lemma [T] and [2j The fact that stationary solution is ergodic is obtained 
by the same argument as in Bougerol and Picard (1992a). 

/s-l \ 

Remark 2 Since 7^ (<I>) < logp I fl 4>s-v ) T-^^e Kesten and Spitzer, 1984j, the condition (jg.jp is neces- 
sary and sufficient for the existence of a strictly stationary solution in of (jS.^P . 

4 Existence of the higher-order moments 

In this section, we present a necessary and sufficient conditions for the existence of finite higher order 
moments for a PGARCH process. 

Theorem 4 Let (1^Jj>q be the stationary solution of model l\2.3^ . Assume K2 < +00 where Kk ='■ E {^t'^}- 
1. // 

p(n^{e\(!/,)}) <i (4.1) 



then Hi G L2. 

2. Conversely, if p ^Yl ^ \^4'f-viVf)^^ — 1> then there is no strictly stationary solution {Y_t)i^'£ to model 
^EM such that E {Xf^] < +00. 



Proof. 
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1. We first define the following M*'^— valued stochastic processes 

Snit) = 



if n < 

B{v,)+A{n)S^^iit-l) ifn>0, 



and for all n S Z: A„(t) = — iS„_i(t). It can be easily shown that, for all n > 0, and 

A„(t) are measurable functions of •••) Hence the processes (<S„(i))tg2 and (A„(t))^g2 

are stationary. From the definition of S_^{t) and A„(t), we can verify that 



Anit) 



if n < 

B{ri^) if n = 

A{n)A^_,{t-l) ifn>0 



For ah n G Z, define T^n\t) = ATi*) = Fee {A„(t)A^(t)}. Using the properties of Kronecker prod- 
uct, we obtain v'nHt) = A'^^ {ri^)r^^l^{t - 1) for n > and E {r^n\t)^ = (^(2))" ^ |r[,2)(t - n)} = 

E [B^Hrj,_j}, where ^(2) := {a®^ (^^)) } . Since p {Ai^)) = p ^11^ {ct>TA%)]) < h 

we conclude that (t) converges in and almost surely to some limit G which is the solution 
of equation (|2.3p . This completes the proof. 

2. From (13.21) we obtain 



oo II k—1 I I oo 

k=0 I I j=0 I I k=0 



^ (2) 



where S^^^ := E |:B®^(r/^)| and the conclusion follows. 

The result of the above theorem can be further extended to the higher-order moments. 

Theorem 5 Let iY_^)^^Q he the stationary solution of model \2.3^ . Assume that H2{r-i) < +oo where 
r > 2. 



1. //p(n^i?{'/'f-.(!Zi)}) < 1 thenY,el.r 



2. Conversely, if p(^^Y\ -E' J ^ 1; then there is no strictly stationary solution (Y_t)t^'j^ to model 
([13) such thatE{Y_f''] < +oo. 

Proof. Define 5:„(t) and A„(i) as in TheoremHand let T'>n\t) = A®'(t). Since ri''^(t) = ^®''(r/^)r[2i(i-l) 
for n > 0, the proof is similar to that for the Theorem H] and thus we omit the details. 
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5 Covariance structure 



To get the covariance structure of the squared PGARCH process, we first assume that the Condition (|4.ip 
holds, this imphes that (|2.2p has a unique (in sense) PC solution. Taking expectation on both sides of 
(12. 2p and using the notation /i, (f ) = E \ y'^j^ >, A; = 1, 2 gives 

ttiiv) = (t)y[i^{v - I) + B^, V = 1, ...,s. (5.1) 

Iterating (jS.ip s— times and requiring ^^(0) = /U^(s), we obtain 

/ s-i \ -1 s-l (j-i ^ 

= I{d) - n 'i^^-'" ^'^-^ \ ^-r 

\ v=0 / j=0 {v=0 ) 



Further manipulations in (jS.ip provide the desired seasonal mean 

j=0 u=o J U=o J 

The seasonal variance can be obtained as follow 

= #/f2(t^-l)+^, 

where ip^ = + (i? {^.(r/^) ® ^,(^^) + 5,(r/^) ./..(r?^)}) /i^(t; - 1) with ^l'^ := E[ct>T{%)} and 

Sl') =^{Bf (^^)}. Thus 

/ s-i \ ~i s-l ci-i 

*,(2) \ J TT ^(2) 



\ v=0 J j=0 U>=o 



j=o u=o ) U=o J 

Now, noting that for any h > 0, we have 

h-1 p-1 ^ r/i-i ^ 

= Y.\ll ('^^-^ ® -^W) f ^v-k®li-y{v -h-k) + \W {4>^_i 1(d)) \ H^iv - h). 
k=0 U=o J U=o J 

In the above equations, cj)^, fj,^{v) and tJ-^i^) are interpreted periodically in v with period s. 



11 



Example 2 The simplest, but often very useful PGARCH process is the PGARCH{1, 1) process. Further 

s s 

computations show that (13. jp and i \4-i\) reduce to Yl ^i(^) < 1 o,nd Yl G2{v) < 1 respectively where 

v=l v=l 

6i{v) = ai{v) + Pi{v) and where 62{v) = K20.i{v) + Pi{v) + 2ai{v)Pi(v) which we henceforth assume. 
However, it is worth noting that the existence of explosive regimes ( i.e. 6i{v) > 1) does not preclude the 
periodic stationarity. The seasonal moments are 
/ s \ -1 s-1 rj-1 ^ 

V v=l / j=0 U=0 J 

j=0 U=o J u=o J 

V v=l J j=0 U=0 J 

^2{v) = E 1 n ^2(^ - ^ {«o(^ - i) + 2«o(^^ - j)ei{v - j)fi,{v - j - 1)} + m 92{V - i) \ fi^is). 

j=Q [i=0 ) [i=0 ) 



The covariance structure of PGARCH (1,1) process can be obtained as 



ao{v)ni{v - 1) + 9i{v)fi2{v -1), h = l 
ao{v)fi,{v - 1) + 9i{v)j,_,{h -l),h> 2. 



When compared to the covariance function of second order GARCH{1, 1) process, the above formulas are 
quite complex. On can verify that these expressions reduce to the classical GARGH{1, 1) forms when 
the PGARCH (1, 1) parameters are constant in v. In general, calculation in PGARCH process are very 
tedious. 



6 Geometric ergodicity and strong mixing 

The basic tools presented here are drawn from the monograph by Doukhan (1994). In this section we 
analyze the statistical properties of PGARCH process, such as the geometric ergodicity and the strong 
mixing. These concepts are fundamental in central limit theorem and law of large numbers which can be 
employed to derive asymptotic normality, consistency of maximum likelihood style estimators and inference 
with the model. 

Definition 1 Let (X^)^^q be a discrete Markov chain taking values in M*^ k > 1 with time homogeneous 
t-step transition probabilities i.e., P* (x, A) = P (X^ € ^|:^o ~ — ) where x € M'^, A ^B^k and where B^k is 
a Borel a— field on with P^ (x, A) = P (x, yl) . Let vr be the invariant probability measure on (M'^, , 
i.e.,y A eB^k: tt{A) = [ P{x, A)TT{dx). 
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a. The chain {2[.t)t>o ■^'^^^ ^'^ geometrically ergodic, if there exists some < r < 1 such that 

VxGM^ \\P\x,.)-^{.)\\y = o{r^) (6.1) 

where \\.\\y is the total variation norm. Furthermore, if the chain (XJ^^^ is started with an initial 
distribution tt, the process is strictly stationary. 

b. The (5— mixing coefficients are defined by 

Pxik) = e[ sup \P{B\a (X„ t < 0)) - PiB) 1 1 = ^ sup ^ \P{A, n B,) - P{A,)P{B,)\ 

where the supremum of the last equality is taken over all pairs of partitions {Ai, Aj} and {Bi, Bj} 
of Q such that Ai ^ a (Xf,t < 0) for all i and Bj £ a {X_t,t > k) for all j. The process (Xf)^^^ is 
called P -mixing if lim Pxi^) — 0- 

This powerful result helps us to develop asymptotic results, since convergence in distribution is ensured for 
all measurable functions of the chain. Perhaps one of the most well-known criterion used in establishing the 
geometric ergodicity of a Markov chain is the drift condition developed in Tweedie (1975) and employed 
for the analysis of stochastic stability. 

Definition 2 A drift function g : M.^ —>■ [l,oo[ satisfies a 1-step geometrical drift criterion (relative to a 
Markov chain) if there exists a compact K CM.^ and a positive constants b and < A < 1, 



E{g{X,)\X,_,=x}< 



A5(x) ifx^K 
h if x G K 

where g is interpreted as a generalized energy function and the compact set K as the centre of attraction. 



Remark 3 One consequence of the geometric ergodicity is that the stationary Markov chain {2Lt)t€Z 1^' 
mixing, and hence strongly mixing, with geometric rate. Indeed, Davidov (1973) showed that for an erqodic 
M...„„ a,,, .... ,r.^m m.as.,. ^ / ||P' ,) - M*). 

13 ^{k) = o{p^) i/dOD holds. 

In what follows, we shall assume, without loss of generality, that p = q, otherwise zeros can be filled. To 
derive the geometric ergodicity and /3— mixing results, we must assume the following conditions for the 
sequence {ri,] 

V-*/ t€Z 

[B.O] The i.i.d. sequence (j]^ has a probability distribution function absolutely continuous with 
respect to the Lebesgue measure and such that the density function take positive values almost surely on 
its support E CW. 

We are now ready to state the main result of this section which establishes statistical properties of 
PGARCH processes, in particular, we focus on the /3— mixing property with exponential decay. 
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Theorem 6 Under Conditions (B.O) and l\3.1^ and assume that 

1. p(n^«) < 1 



\v=l 



2. There are < r < 1 and p < 1 such that E 



then the process {Y_^)^^^ defined by is geometrically ergodic. Moreover, if initialized from its invariant 
measure, (Hj^g^ strictly stationary and (3— mixing with exponential decay. 



Remark 4 The condition E 
E jlog A{ri^) I < anrf I 
has derivative i?'(0) = |lo: 



I < 1 for some < r < 1 in some neighborhood of zero is satisfied if 
} < oo /or some e > 0. Indeed, the function d{v) = E 
I < 0, hence ^{v) decreases in small neighborhood of zero, and since 
T?(0) = \ it follows that i?(r) < 1 for < r < 1. On the other hand £^ | A{r}^ | < 1 implies that 
E |log A{r^^ I < 6t/ Jensen's inequality. 



Proof. The proof is based on Carrasco and Chen (2002, Theorem 1). Assumption B.O ensures that 
Carrasco and Chen's (2002 p. 20) condition on the innovation term holds. We must then prove that 
A {v^ aiid i^t) their assumptions {Aq), {Ai), (Ag) and (A3) 

[A.O] Both A {jl^ and i^t) ^^'^ polynomial functions of r/^, therefore the measurability condition is 
trivially satisfied. 

[A.l] We must show that p {A{Q)) < 1. It turns out however that the nonzero eigenvalues of A{0) are the 
nonzero eigenvalues of $(0). It is not hard to see that then p ($(0)) = p Horn and Johnson, 

1999, p. 68). 



\v=l 



k-1 



[A*. 2] ThatJ] n^U-J^^U 
fc>0 I j=0 



converges almost surely to some constant immediately follows 



from Lemma [TJ It remains to show that < 



fc-i 

j=0 



rj^ .] ^ X converges almost surely to zero for all X G 



fe-i 



as /c ^ cxD. As X is nonrandom, it suffices to show that < A 



7=0 



converges almost surely to zero 



A;-l 



as A; ^ 00. Since E < A 



j=0 



a'', then the condition (13. ip implies that < +cx3 and hence 



fe>0 



fc-i 



lim A^ = 0. Combining this result with the fact that lim < T f (r/ . ) > > yields that almost surely 



k—^oo 



j=0 



14 



fe-1 



0. 



fc— »oo I V~* 1 

\j=0 

[A*. 3] Let z_ € M'*'^ and the Lyapunov function g : 



[1, oo[ by g{z) = + 1, we have 
E{9{Y,)\Y,_,=y} = E{\\Airi^)Y,_, + Biri^)f\Yt_, = y] + l 

mo) 



< E 

< p \\y\\' + 1^ 



y +E 



+ 1 



where k is a positive constant such that E Bij]^ 
M 



< K — 1. Choose A > with 1 — \ > p and set 



and consider the compact K = £ : ||?/|| < M^- 1. For all y ^ K, have (1 — A — p) 

(6.2) 



||y|| > K and therefore 

E {g{Y^)\Y^_^ = y]<P ||y|r + {i-\-p) M < (l - A) g{y). 
When y £ K, we have 



pK 



E{9{Y,)\Y,_,=y}<p\\y\\'^ + ^,<- 

— "-" 1 — A — yO 

Combining (|6.2p and (|6.3p we obtain 



+ K := 6 



(6.3) 



E {9{Y_^)\Y_,_^ = y} < (1 - A) g{y) + hlK^y) 
where \k denotes the indicator function of the K. 

7 Conclusion 

The paper partially extends \? structures of the usual GARCH model to periodic ones which allow the 
volatility of time series to have different dynamics according to the model parameters which switches 
between s-regimes. Our study is based (in a multivariate framework) on a generalized autoregressive rep- 
resentation which we are preferred for a periodic ARM A (PARMA) representation. The main advantage 
of the approach is that, besides its simplicity, it preserve the mathematically tractable GARCH structure 
when s = 1. A thorough examination of the structures of the series and its powers revealed that, under 
appropriate moment conditions, these structures were those of periodic AR (PAR) processes. Beside the 
conditions ensuring the existence and uniqueness of strictly stationary solution of PGARGH , we have 
also gave sufficient conditions for the PGARGH processes to belong to L^, p > 1. Some fundamental 
probabilistic properties such as the /3-mixing and the geometric ergodicity with exponential decay have 
been studied. 
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